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I.  INTRODUCTION 

The  parabolic  equation  method  for  predicting  underwater  acoustic  propagation 
was  introduced  to  the  underwater  acoustics  community  by  Tappert  (1977).  The  parabolic 
equation  (PE)  is  a  one-way  wave  equation  approximation  to  the  full  Helmholtz  wave 
equation  that  defines  the  continuous-wave  (CW),  single  frequency  acoustic  field  in 
general  media.  Because  the  PE  is  still  a  full-wave  equation,  it  includes  all  of  the  full- 
wave  effects  (e.g.,  diffraction)  that  are  neglected  in  more  rudimentary  ray-based  models. 
It  is  also  well  suited  to  handle  range-dependent  media  as  it  naturally  computes  the 
coupling  between  vertical  normal  modes  of  the  waveguide.  Additionally,  the  one-way 
formulation  of  the  PE  allows  for  implementation  with  highly  efficient  numerical 
algorithms.  The  most  common  methods  are  the  split-step  Fourier  (PE/SSF)  method 
(Hardin  and  Tappert,  1973),  the  implicit  finite  difference  (IFD-PE)  method  (Lee,  et  al., 
1981),  and  the  finite  element  (FEPE)  method  (Collins,  1990). 

Of  these  methods,  the  split-step  Fourier  method  was  adopted  in  developing  the 
University  of  Miami  Parabolic  Equation  (UMPE)  Model.  The  UMPE  Model  was 
implemented  by  Smith  and  Tappert  (1993).  Since  Prof.  Smith  moved  to  Monterey  in 
1995,  the  UMPE  Model  was  upgraded  to  incorporate  some  improved  numerical 
techniques  and  rewritten.  This  new  version  was  then  named  the  Monterey-Miami 
Parabolic  Equation  (MMPE)  Model. 

The  MMPE  Model  was  written  in  Fortran  and  implemented  using  a  command  line 
structure  available  in  either  a  Unix  or  DOS  command  window.  A  Fortran  compiler  is 
required  to  create  the  executable  program.  Its  output  file  is  a  binary  data  file  composed 
of  the  complex  values  of  the  acoustic  field  function  for  various  grid  points  in  range, 
depth,  and  azimuth  for  multiple  frequencies,  in  general.  From  this  data,  acoustic  pressure 
and  transmission  loss  can  be  extracted  for  a  variety  of  simulated  receiver  locations  within 
the  computational  grid.  This  post-processing  of  the  binary  data  file  and  the  subsequent 
graphical  manipulation  of  the  results  are  performed  using  Matlab. 
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Many  underwater  acoustic  models  are  written  in  Fortran  language.  The  compiled 
executable  files  run  and  output  data  files  of  the  acoustic  prediction  are  produced.  This  is 
then  processed  using  a  graphics  package  to  display  the  results  of  interest.  Such  command 
line  interface  programs  have  little  advantage  in  Windows  operating  systems.  First  of  all, 
users  have  to  handle  the  programs  separately.  Most  users  would  probably  prefer  using  an 
integrated  software  package  because  it  is  easier  to  handle.  Second,  the  program's 
command  line  interface  and  the  Matlab  command  window  are  text  based.  Users  are 
required  to  type  text  in  the  command  line.  It  is  not  a  user-friendly  computing 
environment.  Third,  in  standard  Fortran  the  size  of  data  arrays  must  be  declared  before 
they  are  used.  Developers  must  always  consider  the  maximum  array  size  that  may  be 
needed.  This  can  sometimes  make  the  programs  waste  unnecessary  memory  depending 
on  some  input  data,  or  users  might  want  to  use  more  input  data  than  the  maximum 
declared.  Fourth,  the  Fortran  program  is  a  stand-alone  program.  It  is  implemented  in  an 
independent  computing  environment.  It  is  not  always  easy  to  combine  with  other 
software  or  systems,  and  it  requires  a  Fortran  compiler  (not  standard  on  all  computer 
systems). 

To  address  these  issues,  the  MMPE  model  must  be  re-programmed.  It  should  be 
an  object-oriented  program,  and  the  propagation  model  and  post-processing  graphical 
program  should  be  integrated  as  one.  It  should  have  a  new  graphic  user  interface  (GUI) 
and  should  be  convertible  to  any  kind  of  system.  For  these  reasons,  the  choice  for  the 
programming  language  is  Java.  Java  is  an  object-oriented  language  and  network 
language.  Java  can  dynamically  create  the  data  array  size  because  it  treats  a  run  time 
object  system.  It  is  also  possible  to  be  multi-threaded,  which  will  allow  it  to  run  several 
tasks  in  parallel  within  the  same  program.  Lately,  the  Java  language  has  become  very 
popular  among  many  programmers  since  it  was  announced  in  1995.  It  is  open-source  and 
freely  served  by  Sun  Micro  Systems  Corporation.  More  and  more  class  libraries  have 
been  and  continue  to  be  developed. 

The  main  goal  of  this  thesis  is  to  develop  a  user-friendly,  efficient,  and  accurate 
full-wave   acoustic   propagation    model    with   enhanced   visualization.      This    will    be 


accomplished  by  integrating  the  MMPE  Model  and  graphic  functions  into  a  single 
program  with  a  new  graphic  user  interface  (GUI).  The  integrated  software  will  make  it 
easier  for  users,  such  as  the  war  fighter  in  an  operational  situation,  to  assess  the  acoustic 
environment  of  the  battle  space. 

Chapter  II  provides  a  summary  of  the  theory  of  the  MMPE  model  and  its 
implementation  and  graphic  manipulation.  It  derives  the  general  form  of  the  parabolic 
equation  to  be  used  and  describes  the  algorithm  used  in  the  MMPE  model.  The  post- 
processing and  graphic  manipulation  routines  are  also  defined.  In  Chapter  in,  the  design 
of  the  Java-based  WinMMPE  and  its  functions  are  presented.  The  details  of  the  GUI  are 
also  presented.  The  model  implementation,  post-processing,  and  graphic  manipulation  are 
described.  In  Chapter  IV,  the  output  results  of  the  WinMMPE  are  described.  Sample 
outputs  of  WinMMPE  are  compared  with  data  from  the  MMPE  Model.  This  chapter  also 
presents  a  discussion  of  the  efficiency  in  terms  of  main  memory  used  for  implementation 
and  running  time.  Finally,  Chapter  V  provides  general  conclusions  and  recommendations 
for  further  work.  Suggestions  for  methods  of  improving  efficiency  and  the  addition  of 
new  modules  are  provided. 
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II.  MONTEREY-MIAMI  PARABOLIC  EQUATION  (MMPE)  MODEL 

A.  INTRODUCTION 

This  chapter  provides  a  general  description  of  the  Monterey-Miami  Parabolic 
Equation  (MMPE)  Model,  its  algorithm,  and  its  implementation.  All  contents  in  this 
chapter  are  summarized  from  the  report  of  the  University  of  Miami  Parabolic  Equation 
(UMPE)  Model  by  Smith  and  Tappert  (1993)  and  the  recent  article  of  Smith  (2000)  that 
analyzed  the  accuracy  of  the  MMPE  Model.  Accordingly,  the  nomenclature  used  in  this 
chapter  is  the  same  as  in  those  reports.  For  more  details,  the  reader  should  refer  to  those 
papers. 

B.  THE  GENERAL  THEORY 

We  begin  by  considering  a  general  time-harmonic,  continuous-wave  (CW) 
acoustic  field  at  angular  frequency  co  in  a  cylindrical  coordinate  system  (r,  z,  (p), 

P(r,z,(p,COt)=  p(r,z,<p)e-,a"    .  (2.1) 

Combining  Eq.  (2.1)  with  the  inhomogeneous  wave  equation  for  a  point  source  (Jensen, 
et  al,  2000), 


pV 


(  1   --"!        1    d2P 


—  VP 


c2    dt 


=  -And  (x  -3c,)  ,  (2.2) 


leads  to  the  Helmholtz  equation  for  the  CW  acoustic  field, 
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In  this  equation,  k0  is  the  reference  wavenumber  defined  as 
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c0  is  the  reference  sound  speed,  n(r,z,#>)  is  the  acoustic  index  of  reference  defined  as 

n(r,z,ff)=      C°        ,  (2.5) 

c(r,z,p) 

and  c(r,z,<p)  is  the  acoustic  sound  speed  that  varies  throughout  space,  in  general.  The 
source  function  is  defined  as  a  point  source  with  reference  source  level  P0  at  a  reference 
distance  of  R0  =  lm,  and 


S{x)  =  ^-S(z-zs)S(r),  (2.6) 

2m 


where  r  =  0  ,  z  =    z  s  refer  to  the  position  of  the  source.    Solutions  to  Eq.  (2.1)  may 

be  computed  over  a  bandwidth  of  frequencies,  and  the  results  can  then  be  Fourier 
synthesized  to  produce  predictions  of  the  arrival  time  structure  due  to  broadband  pulse 
propagation. 

In  most  ocean  environments,  the  influence  of  the  azimuthal  coupling  term  is 
negligible  and  may  be  ignored.  To  further  simplify  Eq.  (2.3),  we  shall  also  neglect  the 
source  term  initially.  We  may  then  define 
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p(r,z)  =  -j=u(r,z)  ,  (2.7) 


and  substitution  into  Eq.  (2.3)  yields  the  free-field  wave  equation 

u  =  0  .  (2.8) 
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Next,  we  introduce  the  operators 
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In  the  far  field  (k0r  » 1 ),  Eq.  (2.8)  then  becomes 
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(2.11) 


Expanding  Eq.  (2.1 1)  into  incoming  and  outgoing  wave  equations,  we  obtain 
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for  the  outgoing  wave  when  the  incoming  wave  is  negligible.   Now  u  can  be  defined  in 
terms  of  the  PE  field  function,  if/(r,  z) ,  as 


w  =  ^(r,z)<?'v  . 

Substituting  this  into  Eq.  (2.12),  we  finally  obtain  the  PE  for  the  field  function, 

-^-=  -ik0y  +  ik0Qw  =  -ik0Hw  , 


(2.13) 


(2.14) 


where 


ff„=i-a. 


(2.15) 


is  an  operator  having  Hamiltonian  form.    This  operator  defines  the  evolution  of  the  PE 
field  function  in  range. 


C.  THE  SPLIT-STEP  FOURIER  (SSF)  ALGORITHM 

We  introduce  the  PE  field  function,  \f/,  by  formally  defining  the  relationship  with 
the  acoustic  pressure  field  as 

P(r,z)  =  P0J^Q-olJ2il/(r,z)e,k°r  ,  (2.16) 

V   r 

normalized  such  that  at  r  =  R^,  \y/\  =  1  and  \p\  =  P0 .  The  computational  evolution  of  the 
field  is  now  accomplished  via  a  marching  algorithm  of  the  form 

y/(r  +  Ar)  =  <$>{r)y/{r)  .  (2.17) 

From  Eq.  (2.14),  the  propagator  3>(r)  can  be  shown  to  have  the  approximate  form 

$(r)  «  e"  Wr)Ar  (2.18) 

where 

H^-^-fdr'Hjr').  (2.19) 


If  the  range  step  Ar  is  small  enough,  the  operator  may  be  assumed  to  be  relatively 
constant  over  a  range  step.  Evaluating  the  operator  at  mid-step,  for  instance,  gives 


Hop(r)  =  Hop 


(       1, 
r  +  —  Ar 
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\ 


(2.20) 


From  Eq.  (2.10),  it  is  observed  that  the  operator  Qop,  and  therefore  Hop,  are 
pseudo-differential  operators  due  to  the  square  root.  An  approximation  must  then  be 
made  for  numerical  implementation.  In  order  to  implement  the  resulting  parabolic 
equation  using  the  split-step  Fourier  method  (Hardin  and  Tappert,  1973),  it  is  necessary 
to  separate  the  index  of  refraction  term  from  the  vertical  derivative  term,  i.e., 

Hop=Top+Uop  .  (2.21) 


It  has  been  shown  (Jensen,  et  al.,  2000)  that  a  solution  which  is  second  order  accurate  in 
range  step  size,  A/%  is  obtained  when  the  propagator  function  is  defined  as 

0(r)  =  e-'k°Tu^r+Ar)e-*o*rroPe-*<>Yu°>(r)  _  (2.22) 

Using  the  wide-angle  operator  splitting  introduced  by  Thomson  and  Chapman  (1983),  we 
define 
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and 

Uop=-(n-l).  (2.24) 

Because  U     is  a  scalar  operator,  it  forms  a  diagonal  matrix  which  can  simply  be 
multiplied  by  the  PE  field  function  at  each  range  step.    However,  Top  is  a  differential 

operator,  is  not  diagonal,  and  couples  the  solution  between  depth  points.   The  split-step 
Fourier  algorithm  overcomes  this  by  recognizing  that  the  corresponding  operator  in  the 

vertical  wavenumber  domain,    Top ,  is  diagonal  and  can  be  treated  using  a  simple 

multiplication. 

In  the  PE/SSF  implementation,  generally,  a  multiplication  of  the  z-space  operator 

Ar 
-ikQ—Uop(r) 

e      -  is  defined  at  the  initial  range-step.    The  result  is  then  transformed  to  the 

vertical  wavenumber,   kz,  domain.     After  a  multiplication  by  the   kz -space  operator 

e  '  °  r  °p ,  the  solution  is  transformed  back  to  the  z-domain  where  a  multiplication  of  the  z- 

-ikQ—Uop{r+Ar) 

space  operator  e      2  is  defined  at  the  end  of  the  range  step.     A  fast  Fourier 

transform  (EFT)  algorithm  is  used  to  accomplish  the  transformation  which  assumes  a 
convention 


¥(z)  =  FFT{v{kzj) 


(2.25) 


The  representation  of  the  PE/SSF  implementation  is  then 


i/f(r  +  Ar,z)  =e 


-«*0^(r+Ar,z) 

2  x  Ft  T 


e-ik0*rTop(kz)  jFFT 


-iko^Uop(r,z) 

e      2  y/(r,z) 


(2.26) 


where,  in  k.  -space, 


T  =  \- 


1- 


fkA 


i1/2 


(2.27) 


D.  THE  BOUNDARY  CONDITIONS 

The  ocean  surface  in  the  MMPE  model  is  defined  by  a  Dirichlet  boundary 
condition  consistent  with  a  pressure  release  boundary, 

V(z  =  0)  =  0  .  (2.28) 

For  this  boundary  condition,  an  image  ocean  method  is  used  whereby  an  exact  copy  of 
the  environment  is  defined  for  negative  depths  and  the  field  function  for  these  depths  is 
defined  by 


y/{-z)  =  -y/(z) 


(2.29) 


Accordingly,  the  data  array  is  twice  as  large  as  that  needed  to  describe  the  real  acoustic 
field. 

In  the  MMPE  model,  shear  wave  propagation  in  the  bottom  is  treated  by  an 
equivalent  fluid  technique.  Thus,  the  propagation  of  shear  waves  are  not  computed  but 
the  effect  of  shear  on  the  boundary  reflection  coefficient  is  included.  Specifically,  the 
equivalent  fluid  approximation  of  Tindle  and  Zhang  (1992)  is  employed  at  each  bottom 
interface. 
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It  is  assumed  that  the  interface  between  the  bottom  of  the  water  column  and  the 
top  of  the  basement,  or  sediment,  layer  is  characterized  by  a  discontinuity  due  to  a  sharp 
contrast  in  sound  speed.  Specifically,  the  sound  speed  at  the  interface  is  defined  by 

c(z)  =  cw(z)  +  [cb(z)-cw(z)]H(z-zb)  ,  (2.30) 

where   cw   and   cb   are  the  sound  speed  of  water  and  bottom,  respectively,  and  the 
Heaviside  step  function  is  defined  by 


H(g)  = 


[0,    (<Q 
1/2,    C  =  0 
l,    £>0 


(2.31) 


where  £  s  z  -  zb  ■  In  order  to  implement  within  the  numerical  scheme,  it  is  necessary  to 

approximate  this  by  a  smooth,  continuous  function.    For  that  purpose,  this  function  is 
approximated  by  a  hyperbolic  tangent  function, 


*<o=| 


1  +  tanh 


2L 


V 


J-i 


(2.32) 


or,  equivalently, 


tf(0  =  (i 


■*)  ' 


,SILc 


(2.33) 


where  Lc  is  the  sound  speed  mixing  length. 

In  order  to  create  the  most  realistic  interface  condition  for  reflections  from  a 
sound  speed  discontinuity,  Smith  (2000)  showed  empirically  that  the  most  accurate 
results  were  obtained  by  defining  Lc  as  a  fraction  of  a  wavelength,  specifically  AJ\Q.  For 
sampling  purposes,  the  default  minimum  value  for  Lc  is  set  at   Lc  miD  =  Az .     Thus, 

Lc  =  max(Az,/l/lO),    so    the    most    accurate    results    are    obtained    efficiently    when 
Az<i/10. 
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We  must  now  consider  the  effect  of  the  bottom  density  contrast,  something  which 
was  neglected  in  the  definition  of  the  wave  equation  given  in  Eq.  (2.2).  A  proper 
treatment  including  density  produces  the  same  expressions  given  above  if  we  replace  the 
index  of  refraction  n  by  an  "effective"  index  of  refraction, 


n     =  n~  +■ 


1 


2k 


1  ^f  1 


p 


p 


(2.34) 


Then,  the  pressure  field  p(r,  z)  is  defined  by 


rtr.d-PoJ—  Q;vV(r,zyv  , 

V  Por 


(2.35) 


where    p0    is  a  reference  density  defined  by    p0=l.O   g/cm" 
propagator  function  in  Eq.  (2.24)  must  now  be  defined  as 

Uop(z)  =  Ul(z)  +  U2(z)  , 


The  environmental 


(2.36) 


where  Ux{z)  is  the  same  as  the  previously  defined  environmental  potential  function  in 
Eq.  (2.24)  and  U2(z)  is  a  new  operator  that  includes  the  effect  of  the  density 
discontinuity. 


If  we  assume  that  p  =  p{z)  only,  the  density  is  defined  by 
P(z)  =  pw+(pb-pw)H{z-zb)  , 


(2.37) 


where  pw  and  pb  are  the  density  of  water  and  bottom,  respectively.  As  before,  we  must 
replace  the  Heaviside  step  function  by  a  smooth,  continuous  function  scaled  by  a  density 
mixing  length,  Lp  .  Within  the  MMPE  Model,  the  function  £/2  is  approximated  by 


A:n  dz 


(2.38) 


where 
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£  = 


i+(pJpJ/2 


(2.39) 


For   the    density   mixing   function,    a   cubic    spline    over   the    finite   interval 
L    <  £  <  Lp  is  defined.  Then  H(£)  is  defined  by 


H(C)  = 


0 


2     Lf     3 


1  +  — 


<r 


3 


f 


v    p  J 


1- 


V 


L, 

<c<- 

2 
L, 


-Lp<-         -^r 


~^C 


2  ; 


OK 


(2.40) 


The  first  derivative  of  this  function  is 


H  (O  =  S(C)  = 


_2 


C 


1-2 


0 
1  +  -2- 


vL^ 


1-X 
0 


_L  <r<-hiL 

'         2  "S"      2 
'         2   ~  4  "    ' 


(2.41) 


where  S(£)   is  the  smooth  approximation  to  the  Dirac  delta  function.     The  second 
derivative  is  then 
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H  (0  =  S(0  = 


_4_ 
L\ 


0 

(       C 

1  +  -2- 

Lp 


P 
_4 


f 


\ 


y 


V    V 


1- 


<r 


V 


0 


C^-L 


p 
L, 


,    -L   <<T<--f 


2  2 


<T^ 


(2.42) 


Equations  (2.38)  and  (2.42)  provide  the  necessary  expressions  for  computing  the  density 
potential  function  with  a  cubic  spline  polynomial  smoothing  function. 

In  order  to  create  the  most  realistic  interface  condition  for  reflections  from  a 
density  discontinuity,  Smith  (2000)  showed  empirically  that  the  most  accurate  results 
were  obtained  by  defining  L    ~  2A  .  For  sampling  purposes,  the  default  minimum  value 

for  Lp  is  set  at  LpmiII  =  5Az  .  Thus,  Lp  =  max(5Az,2/l),  so  the  most  accurate  results  are 

obtained  efficiently  when  Az  <  2A/5. 

E.  THE  SOURCE  FUNCTION 


Near  the  source,  we  may  assume  that  the  effect  of  density  is  negligible  and  the 
operator  in  Eq.  (2.16)  is  near  unity.  The  acoustic  pressure  near  the  source  is  then  defined 
simply  by 


\R 


0  ,„„iV 


p  =  P    -a-¥e 
r 


(2.43) 


where  R0  is  a  reference  range  defined  by  lm  and  P0  is  the  source  pressure  at  R0.   The 
source  level,  SL,  becomes 


5L  =  201og 


dB  re  Pr,RQ  , 


(2.44) 
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where  Pr  is  a  reference  pressure,  typically  defined  by  Pr  =  1  juPa 
From  Eq.  (2.43),  y/(r  =  0,  z)  can  be  expressed  as 


1       r 
y/(r  =  0,  z)  =  lim—    — p(r,  z)  .  (2.45) 


o  V     o 


In  the  vicinity  of  a  point  source,  the  pressure  field  takes  the  form  of  a  spherical  Green's 
function,  i.e. 

p=—e**     ,    R  =  Jr2+z2     ,    a  =  P0R0,  (2.46) 

AnR 

p 

where  a  is  defined  by  requiring    \p\  =—  at  /?  =  P0.    If  we  then  define  the  PE  field 

An 

starting  condition  for  a  point  source  at  depth  zs  by 

y/{r  =  0,z)  =  a5{z-zs)  ,  (2.47) 

then  it  can  be  shown  that  the  proper  normalization  is 


a=&-  •    {2A8) 


lo 


To  represent  this  field  in  the  vertical  wavenumber,  &z-domain,  we  perform  a  Fourier 
transformation  of  Eq.  (2.47)  to  produce 

fi(r  =  0,k)  =  -2iasm(kzzs)  ,  (2.49) 

which  also  includes  the  influence  of  the  image  source. 

For  the  wide  angle  approximation  in  the  vertical  wavenumber  domain,  we  need  to 
add  two  factors  in  Eq.  (2.49).  First  of  all,  to  produce  the  correct  solution  in  the  far-field, 

f  i  2  V1/4 


the  factor 


Kf. 


,    (\kz\  <k0)  needs  to  be  added  to  the  wide  angle  source  (Thomson 
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and  Bohun,  1988).     Secondly,  in  order  to  take  advantage  of  the  MMPE  depth  grid 
convention,  defined  by 


z„  = 


( 


1 


n 

N  -n  +  — 
2J 


i   N 
Az  ,    n  =  1,  — 

2 


N         1       A7 

Az  ,     n  =  —  +  l,N 
2 


(2.50) 


Az 


we  need  to  add  a  phase  term  in  the  wavenumber  domain  of  e     2  .   Finally  then,  the  PE 
starting  field  in  the  ^-domain  for  the  wide  angle  source  is  given  by 


i/(r  =  0,*.)  =  -2i\ 


iR, 


27ikr 


sin(k.zs) 


2  N 


'i 


-1/4 


Az 


(2.51) 


F.  IMPLEMENTATION 

There  are  two  aspects  of  the  implementation  of  the  MMPE  Model.  The  first  is  the 
propagation  model  implementation  that  produces  a  binary  output  file.  The  second  aspect 
is  then  the  post-processing  of  this  binary  data  file  and  the  graphical  representation  of  the 
results.  This  is  treated  using  Matlab  routines.  The  results  of  the  output  PE  solution  will 
be  discussed  in  Chapter  IV.  In  this  section,  the  run-time  process  and  the  output  structure 
is  discussed. 

1.  Model  Implementation 

The  MMPE  Model  can  be  run  using  a  DOS  command  line  interface  as  illustrated 
in  Figure  2.1  (a).  Prior  to  execution,  several  ASCII  input  files  must  be  defined.  The 
primary  input  data  file  is  named  "pefiles.inp".  Within  this  file,  all  other  filenames  that 
provide  details  of  the  environment  (sound  speed  profiles,  bathymetry,  etc.)  are  defined. 
Of  all  file  names,  only  the  specific  name  "pefiles.inp"  is  required.  Also  within  this  main 
file  are  specifications  for  the  output  data  file  (number  and  range  of  output  depths  and 
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ranges)  as  well  as  optional  declarations  of  computational  parameters  (FFT  size,  range  and 
depth  grid  sizes,  etc.). 

When  run,  the  propagation  model  asks  the  user  to  choose  either  an  accurate  mode 
or  an  efficient  mode  of  calculation  (unless  the  grid  sizes  have  been  explicitly  specified  in 
"pefiles.inp").  In  the  former,  the  depth  grid  is  1/10  of  a  wavelength  and  the  range  grid  is 
a  wavelength.  In  the  latter,  the  depth  grid  is  half  a  wavelength  and  the  range  grid  is  three 
times  a  wavelength.  The  model  calculates  the  PE  field  depending  on  the  user's  choice. 
During  the  calculation,  the  range  step  displays  on  the  screen,  along  with  the  current 
frequency  and  radial  bearing  being  computed,  as  illustrated  in  Fig.  2.1  (b).  The  result  of 
the  calculation  is  the  binary  output  with  filename  defined  in  "pefiles.inp." 

The  structure  of  the  output  binary  file  is  shown  in  Figure  2.2.  For  all  but  the  final 
records  (containing  bathymetry  information),  each  record  length  is  twice  the  number  of 
points  stored  in  depth.  The  factor  of  two  accounts  for  the  real  and  imaginary  parts  of  the 
complex  PE  field  function.  The  first  record  contains  header  information  providing  details 
about  the  remainder  of  the  data  and  the  calculation  as  follows:  the  total  length  of  the  data 
file;  the  reference  sound  speed;  the  number  of  frequencies  computed;  the  center 
frequency;  the  frequency  bandwidth;  the  number  of  range  points  output;  the  minimum 
range  for  output;  the  maximum  range  for  output;  the  range  step  size  used  in  the 
calculation;  the  number  of  depth  points  output;  the  minimum  depth  for  output;  the 
maximum  depth  for  output;  the  number  of  radials;  the  angular  difference  between  radials; 
the  bottom  depth  at  the  source  location;  the  deep  bottom  depth  at  the  source  location;  and 
the  source  depth. 

For  multiple  frequency  and  multiple  radial  calculations,  the  model  first  computes 
the  solution  out  in  range  for  the  initial  frequency  and  radial.  This  begins  by  computing 
the  starting  field,  which  is  the  same  for  all  radials  at  that  frequency.  It  then  computes  the 
solution  for  all  subsequent  radials  at  that  frequency.  After  computing  the  solution  out  in 
range  for  every  radial,  it  then  begins  the  calculations  for  the  subsequent  frequencies.  The 
result  of  the  ordering  is  the  binary  data  file  organization  depicted  in  Figure  2.2. 
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Figure  2.1:  The  MMPE  model  run-time  environment:  (a)  initiation  of  model  run; 
(b)  model  output  while  running. 
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number  of  depths 


frequency  1 


frequency  2 


i 

i                                                            ' 

header 

range  1 
range  2 

riririririri  ....    Starting  field 

radial  1 
radial  2 

riririririri .... 

riririririri  .... 

radial  1 
radial  2 

riririririri .... 

riririririri .... 

r:  real  value 

i:  imaginary  value 


Figure  2.2:  The  structure  of  the  MMPE  output  binary  data  file. 

Note  that  the  data  ordering  places  adjacent  radials  next  to  each  other  rather  than 
adjacent  range  steps  as  might  be  expected  from  the  actual  run-time  ordering.  This  was 
done  to  maintain  consistency  with  other  versions  of  MMPE,  which  incorporate  true  3D 
effects  with  azimuthal  coupling.  In  those  versions,  all  radials  are  computed  at  each  range 
step.  Also  note  that  the  complex  PE  field  function  values  are  stored  as  real  and 
imaginary  data  pairs. 

2.  Post-Processing  and  Graphic  Manipulation 

The  post-processing  implementation  runs  in  the  Matlab  command  window.  The 
initialization  routine  is  called  "peoutl",  which  requests  the  user  to  enter  the  name  of  the 
binary  output  file  of  the  MMPE  Model  to  be  processed.  After  entering  the  file  name,  the 
data  of  the  output  file  is  loaded  to  the  variables  in  "peoutl"  as  shown  in  Figure  2.3  (a). 
The  main  processing  routine,  "peout2",  can  then  be  run  multiple  times  without  re-running 
the  initialization  routine  so  long  as  the  same  file  is  being  processed.     After  typing 
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"peout2",  six  options  for  processing  are  given,  as  shown  in  Figure  2.3  (b).  The  user  then 
chooses  the  option  for  the  desired  processing  and  data  extraction.  For  the  option  chosen, 
all  acoustic  quantities  of  interest  are  computed,  and  a  corresponding  transmission  loss 
field  is  plotted.  The  outputs  will  be  discussed  in  more  detail  in  Chapter  IV. 


For  all  output  options  of  interest  here,  transmission  loss  (TL)  is  defined  by 


TL(r,z,(p)  =  20\og 


10 


\p(r,z,<p)\ 


dB  re  lm , 


(2.52) 


where  Po  is  the  source  pressure  amplitude  measured  at  the  reference  distance  of  R0  =  1  m. 

From  the  definition  of  the  pressure  field  in  terms  of  the  PE  field  function  in  Eq.  (2.16), 
this  may  be  rewritten 

(2.53) 


TL(r,z,<p)  =  201og10 


-=\y/(r,z,<p)\ 
Vr 


=  201og]0|^(r,z,p)|-101ogI0r  dBrelm 

Thus,  the  solution  of  the  PE  field  function  provides  all  the  information  necessary  to 
compute  TL  with  the  exception  of  the  simple  cylindrical  spreading  factor. 
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Figure  2.3:  Post-processing  and  graphic  implementation:  (a)  loading  the  MMPE  output 
file;  (b)  choosing  one  of  the  processing  options. 
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III.  WINDOWS  MMPE  (WINMMPE) 

A.  INTRODUCTION 

As  described  in  the  previous  chapter,  the  MMPE  model  can  be  executed  via  Unix 
shell  or  DOS  command.  Then,  the  post-processing  and  the  graphic  manipulation  are 
implemented  in  MATLAB.  Since  the  final  goal  of  this  thesis  is  integrating  those 
functions  and  providing  a  better  interface  for  users,  a  Windows  version  of  the  MMPE 
satisfying  these  goals  has  been  created.  It  is  called  WinMMPE.  This  chapter  describes 
its  design  and  functions.  Also,  its  implementation  is  described  in  terms  of  the  model  and 
the  graphic  manipulation.  The  WinMMPE  application  is  executed  by  clicking  its  icon 
(Figure  3.1). 


Figure  3.1  -  Icon  of  WinMMPE 


B.  DESIGN  AND  FUNCTIONS 

WinMMPE  is  composed  of  five  frames,  consisting  of  the  main  frame,  model 
control  frame,  input  data  editor  frame,  graphic  control  frame,  and  graphic  output  frame. 
Each  of  the  frames  has  its  own  design  and  functions.  This  section  describes  how  they 
have  been  designed  and  what  kinds  of  functions  they  have. 

1.  Main  Frame 

A  double  click  of  the  left  mouse  button  on  the  icon  makes  the  WinMMPE  display 
the  main  frame  (Figure  3.2).  It  provides  a  full  Windows  graphic  user  interface  for  all 
WinMMPE  functionality.  There  are  five  menus,  which  are  "File"  menu,  "Edit"  menu, 
"View"  menu,  "Window"  menu,  and  "Help"  menu.  Each  menu  has  several  menu  items 
and  all  are  shown  in  Table  3.1.    In  the  "File"  menu,  there  are  five  menu  items.    "New 
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Source"  menu  item  renders  model  control  frame  with  main  input  file,  "default.inp." 
"Open  Source"  menu  item  renders  an  open  file  chooser.  Users  can  choose  a  "inp"  file 
as  a  main  input  file.  If  other  formats  are  chosen  instead  of  ".inp,"  file  errors  will  occur. 
"Save"  and  "Save  As"  menus  are  used  to  save  the  main  input  file  with  the  given  file 
name,  or  to  save  it  with  another  name  if  desired.  "Exit"  menu  item  makes  WinMMPE 
terminate.  The  "Edit"  menu  has  three  menu  items:  "Cut,"  "Copy,"  and  "Paste."  When 
editing  input  files  in  the  "Input  File  Editor"  frame,  these  menu  items  are  supposed  to 
work  for  the  frame.  However,  they  do  not  work  yet  and  are  left  for  future  work. 
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Figure  3.2  -  WinMMPE  Main  Frame 
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In  the  "View"  menu,  there  are  three  menu  items:  "Model  Control"  frame,  "Input 
File  Editor"  frame,  and  "Graphic  Control"  frame.  Whenever  any  of  these  menu  items  is 
selected,  a  file  chooser  is  opened.  For  "Model  Control"  frame  and  "Input  File  Editor" 
frame,  it  shows  ".inp"  file  as  option.  For  "Graphic  Control"  frame,  it  shows  ".bin"  file  as 
option.  The  "Window"  menu  shows  the  opened  frames  as  menu  items.  Users  can  check 
how  many  window  frames  are  on  the  main  frame  and  can  pop  up  the  frame  they  select. 
The  "Help"  menu  has  two  menu  items:  "Web  Help"  and  "About  WinMMPE."  The 
former  lets  users  give  the  information  about  menu  and  each  menu  item  of  WinMMPE. 
The  latter  is  the  general  explanation  of  WinMMPE. 


The  nine  buttons  in  the  tool  bar  below  the  menu  bar  are  seen  in  Figure  3.1.  From 
left  to  right,  they  are:  "New  Source"  button,  "Open  Source"  button,  "Save"  button,  "Cut" 
button,  "Copy"  button,  "Paste"  button,  "Model  Control"  frame  button,  "Input  Data 
Editor"  frame  button,  and  "Graphic  Control"  frame  button.  Each  button  is  able  to  work 
like  the  corresponding  menu  item  previously  described. 

25 


2.  Model  Control  Frame 

This  frame  launches  the  main  input  data  file  and  the  sound  source  input  file, 
which  have  the  same  form  as  the  ones  used  in  the  MMPE  model.  It  comprises  seven 
panels:  "Output  Data  Files"  panel,  "Depth  Info"  panel,  "Range  Info"  panel,  "Basic  Info" 
panel,  "Source  Info"  panel,  mode  selection  panel,  and  buttons  panel  (Figure  3.3).  The 
first  four  panels  load  the  main  input  data  file,  and  the  "Source  Info"  panel  loads  the  sound 
source  input  file.  In  the  "Output  Data  Files"  panel,  the  output  file  name  and  log  file  name 
can  be  given.  The  text  fields  in  the  "Depth  Info"  and  "Range  Info"  panels  decide  the 
depth  and  range  information  for  the  output  data.  The  text  fields  in  the  "Basic  Info"  panel 
provide  parameters  for  running  the  model,  for  example,  FFT  size.  If  the  value  of 
"Number  of  Depth  Points"  is  zero,  depth  grid  size  is  determined  according  to  accurate  or 
efficient  mode.  If  it  is  not  zero,  it  is  assigned  to  the  variable,  deciding  the  FFT  size.  The 
value  of  "Absolute  Max  Depth"  is  the  maximum  depth  of  the  deep  bottom.  If  this  value 
is  zero,  it  is  reset  to  the  two  times  the  maximum  depth  of  the  water/bottom  interface. 

The  radio  buttons  in  the  mode  selection  panel  for  depth  and  range  are  used  in 
deciding  the  grid  size  of  depth  and  range.  The  buttons  panel  has  five  buttons:  "Run" 
button,  "Plot"  button,  "Input  File  Info"  button,  "Clear"  button,  and  "Help..."  button.  The 
"Run"  button  starts  executing  the  model.  After  the  model  computations  are  completed, 
the  output  binary  file  is  created  and  the  "Graphic  Control"  frame  is  automatically  opened. 
The  "Plot"  button  renders  to  open  the  "Graphic  Control"  frame  with  the  output  binary  file 
given  in  the  text  field  of  the  output  file  name.  The  "Input  File  Info"  button  renders  to 
open  the  "Input  File  Editor"  frame  with  the  input  files  recorded  in  the  main  input  file. 
The  "Clear"  button  clears  all  text  fields.  The  "Help..."  button  shows  the  explanation 
about  this  frame. 

3.  Input  Data  Editor  Frame 

This  frame  loads  five  input  data  files,  specifically  sound  speed  profile  (SSP)  data 
file,  bottom  bathymetry  data  file,  bottom  property  data  file,  deep  bottom  bathymetry  data 
file,  and  deep  bottom  property  data  file.    Like  the  main  input  file  and  the  sound  source 
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input  file,  these  files  have  the  same  form  as  the  ones  used  in  the  MMPE  model.  The 
format  of  the  data  files  is  described  in  Smith  and  Tappert  (1993).  The  "Input  Data 
Editor"  frame  has  a  tabbed  pane  with  five  panels  as  shown  in  Figure  3.4.  Each  panel  has 
a  file  name  text  field,  "Save"  button,  "..."  button  which  renders  a  file  chooser  to  open  an 
input  data  file,  and  a  text  area. 
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Figure  3.3  -  Model  Control  Frame 


User  evaluation  found  that  the  "Input  File  Editor"  frame  does  not  work  well  as  an 
ASCII  file  editor.  JTextArea,  a  Java  class  in  the  Java  standard  library,  provides  an 
editing  capability  for  the  input  files.  However  using  this  class  is  not  desirable  because  it 
sets  the  size  of  rows  and  columns  with  two  dimensional  character  arravs.    Therefore, 
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future  work  needs  to  avoid  using  JTextArea  class.    JEditorPane  might  well  be  used 
instead.  The  problems  in  this  frame  are  left  for  future  work. 
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Figure  3.4  -  Input  File  Editor  Frame 

4.  Graphic  Control  Frame 

As  described  in  Chapter  II,  in  the  post-processing  of  MMPE.  there  are  six  kinds  of 
options  for  graphic  manipulation.  In  WinMMPE,  four  options  have  been  worked  so  far: 
"Single  Radial,"  "Single  Depth,"  "Single  Interface"  and  "Source  Field."  The  "Graphic 
Control"  frame  has  a  panel  and  a  tabbed  pane  and  is  shown  in  Figure  3.5.  The  output 
binary  file  is  chosen  in  the  panel.  With  the  "..."  chooser  button  the  output  binary  file  that 
users  want  to  use  can  be  selected  and  loaded.  The  tabbed  pane  has  the  four  panels  with 
each  graphic  option.  The  graphic  output  from  each  graphic  option  and  the  components  of 
each  panel  in  details  are  described  later  in  this  chapter. 
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Figure  3.5  -  Graphic  Control  Frame 

5.  Graphic  Output  Frame 

The  above  three  frames,  "Model  Control"  frame,  "Input  File  Editor"  frame,  and 
"Graphic  Control"  frame  are  displayed  inside  the  main  frame,  while  the  "Graphic 
Output"  frame  is  displayed  independently  on  the  main  frame  and  is  shown  in  Figure  3.6. 
The  title  bar  of  the  frame  shows  the  chosen  graphic  option  from  the  "Graphic  Control" 
frame  and  the  output  file  name.  There  are  four  menus:  "File,"  "Edit,"  "Window,"  and 
"Help."  As  shown  in  Figure  3.6  (a),  the  "File"  menu  has  seven  menu  items:  "New," 
"Open,"  "Save."  "Save  As,"  "Print,"  "Page  Setup,"  and  "Close."  Of  these  menu  items, 
"Print,"  "Page  Setup,"  and  "Close"  do  work  and  the  others  do  not  work.  In  addition, 
"Edit"  and  "Window"  menus  do  not  work,  either.  They  are  left  for  future  work.  The 
"Help"  menu  shows  the  user's  choice  from  the  option  in  the  "Graphic  Control"  frame.  It 
is  shown  in  Figure  3.6  (b). 
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Two  kinds  of  graphic  output  can  be  launched  in  this  frame.  One  is  a  color  map 
image  of  a  transmission  loss  field  (Figure  3.6  (a)).  If  the  user  would  like  to  know  the 
value  of  a  transmission  loss  at  a  point  inside  the  image,  the  user  can  put  the  cursor  at  that 
position  and  then  click.  The  value  is  shown  at  the  bottom  label  of  the  graphic  frame.  The 
bar  of  color  index  with  respect  to  transmission  loss  is  also  given.  Another  type  of  plot  is 
a  single  line  of  the  transmission  loss  (Figure  3.6  (b)).  The  number  and  the  value  of  the 
scale  in  the  x-axis  and  y-axis  of  the  plot  are  decided  by  the  value  depending  on  the 
difference  of  minimum  value  and  maximum  value. 

To  develop  the  algorithm  for  the  color  map  shown  in  Figure  3.6  (a),  a  lot  of 
experimental  trial  and  error  regarding  treatment  of  RGB  values  was  performed.  The  final 
color  method  is  shown  in  Figure  3.7.  The  red,  green,  and  blue  values  of  the  color  map 
start  from  0.5,  0,  and  0,  respectively  at  the  minimum  TL.  They  finish  at  0,  0,  0.5  at  the 
maximum  TL. 
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Figure  3.6  -  Graphic  Output  Frame 
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C.  IMPLEMENTATION 

As  described  in  the  previous  section,  all  of  the  functionality  in  WinMMPE  is 
initiated  by  clicking  a  button  in  the  main  frame.  Model  execution  is  initiated  by  clicking 
the  "Run"  button  in  the  "Model  Control"  frame.  The  graphic  manipulation  is  initiated  by 
the  "Plot"  button  in  the  "Graphic  Control"  frame.  For  graphic  image  generation,  each  of 
four  graphic  options  is  described  in  this  section.  This  section  does  not  describe  what 
kinds  of  data  used,  but  rather  how  to  manipulate  the  graphic  outputs. 

1.  Model  Execution 

As  soon  as  the  model  starts  running,  processing  outputs  appear  in  the  DOS 

command  window  as  shown  in  Figure  3.8.  The  program  then  produces  a  binary  output 

file.   The  file  has  the  same  structure  as  the  one  from  the  MMPE  model.   Ordinarily  the 

binary  serialization  file  in  Java  uses  big-endian  floating-point  representations  as  a  default 

type  in  the  standard  Java  library  supported  by  Sun  Microsystems  Corp.   However,  in  the 
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WinMMPE,  the  little-endian  representation  is  used  for  the  binary  output  file  to  be 
compatible  with  the  Windows  version  of  Matlab  for  post-processing  of  the  MMPE  model 
data.  The  little  endian  Java  class  is  referenced  at  http://mindprod.com/products.html  - 
LEDATASTREAM. 
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Figure  3.8  -  Model  Running  in  the  WinMMPE 
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2.  Post-processing  and  Graphic  Manipulation 

The  "Single  Radial"  post-processing  option  in  the  "Graphic  Control"  frame  is 
shown  in  Figure  3.5.  The  graphic  output  is  produced  by  Eq.  (2.53).  It  has  a  "Radial" 
section,  a  "Frequency"  section,  and  four  buttons:  "Plot,"  "VRML,"  "Clear,"  and  "Help." 
Each  section  has  an  input  text  field.  Whenever  the  output  file  is  loaded,  0.0  degrees  and 
the  center  frequency  are  set  as  defaults  in  the  "Radial"  and  "Frequency"  fields, 
respectively.  The  "Plot"  button  renders  a  "Single  Radial"  graphic  output.  The  title  bar  of 
the  frame  shows  the  phrase,  "Single  Radial,"  and  the  output  file  name.  The  acoustic 
transmission  loss  field  at  a  single  radial  along  the  depth  and  the  range  is  shown  in  Figure 
3.9.  The  single  line  is  the  bottom  bathymetry  as  a  function  of  range.  It  is  possible  to 
show  the  deep  bottom  bathymetry  as  well,  if  it  exists  in  the  data. 
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Figure  3.9  -  Single  Radial  Graphic  Output 
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The  "VRML"  button  is  only  in  the  "Single  Radial"  option.  This  button  renders  to 
create  four  ASCII  data  files  for  showing  a  VRML  scene  graph.  They  are  "point"  for 
point  data,  "color"  for  color  data,  and  "indexPoint"  for  coordlndex  and  colorlndex  in 
VRML.  These  files  are  combined  in  a  VRML  file,  "mmpeVRML.wri."  The  WinMMPE 
then  renders  a  web  browser  set  as  default  on  the  computer.  The  web  browser  must  be 
able  to  launch  the  "VRML"  file  with  the  extension  "wrl."  Then  the  3D  scene  of  the 
acoustic  transmission  loss  field  is  simulated,  as  depicted  in  Figure  3.10.  In  fact,  the  data 
can  be  viewed  with  the  true  aspect  ratio  in  terms  of  real  range  and  real  depth.  Note  that 
the  scene  of  Figure  3.10  is  rninimized  by  1/1000  for  range  and  by  1/100  for  depth. 
Currently  the  3D  scene  does  not  have  much  meaning  since  it  is  simply  a  single  2D 
acoustic  field  in  3D  space.  Future  work  expanding  the  3D  model  is  discussed  in  Chapter 
V. 
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Figure  3.10  -  Acoustic  Field  in  VRML 
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The  GUI  of  the  next  option,  "Single  Depth",  is  shown  in  Figure  3.11.  The 
graphic  output  is  produced  by  Eq.  (2.53).  It  has  a  "Depth"  section,  a  "Frequency" 
section,  and  three  buttons.  When  the  output  file  is  loaded,  the  source  depth  and  center 
frequency  are  set  as  defaults  in  the  "Depth"  and  "Frequency"  fields,  respectively.  The 
"Plot"  button  renders  a  graphic  output  of  "Single  Depth."  Figure  3.12  shows  two  cases 
of  graphic  outputs.  Figure  3.12  (a)  shows  a  multi-radial  acoustic  field  and  Figure  3.12  (b) 
shows  a  single  radial  acoustic  field.  Note  that  if  multiple  radials  exist  in  the  data  file,  the 
plot  automatically  displays  the  TL  field  of  radial  versus  range. 
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Figure  3.1 1  -  The  Option  "Single  Depth"  in  Graphic  Control  Frame 
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Figure  3.12  -  Single  Depth  Graphic  Output:  (a)  Multi-Radial  Acoustic  Field;  (b)  Single- 
Radial  Acoustic  Field 
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The  GUI  of  the  "Single  Interface"  is  shown  in  Figure  3.13.  The  graphic  output  is 
produced  by  Eq.  (2.53).  It  has  a  "Depth"  section,  a  "Frequency"  section,  and  three 
buttons.  Unlike  the  "Depth"  section  in  the  "Single  Depth"  graphic  option,  it  has  two 
radio  buttons  for  choosing  one  of  either  "Water/Bottom  Interface"  or  "Bottom/Basement 
Interface,"  instead  of  an  input  text  field.  When  an  output  file  is  loaded,  "Water/Bottom" 
and  center  frequency  are  set  as  defaults  in  the  "Depth"  and  "Frequency"  fields, 
respectively.  The  "Plot"  button  renders  a  graphic  output  of  "Single  Interface."  Figure 
3.14  (a)  shows  a  multi-radial  acoustic  field  at  "Water/Bottom  Interface"  and  Figure  3.14 
(b)  shows  a  single  radial  acoustic  field  at  "Bottom/Basement  Interface,"  as  shown  in  the 
title  bar  of  each  figure. 
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Figure  3.13  -  The  Option  "Single  Interface"  in  Graphic  Control  Frame 
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Figure  3.14-  Single  Interface  Graphic  Output:  (a)  Multi-Radial  Acoustic  Field;  (b) 
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Figure  3.15  shows  the  GUI  of  the  "Source  Field."  The  graphic  output  is  produced 
by  Eq.  (2.53).  It  has  only  "Plot"  and  "Help..."  buttons.  The  "Plot"  button  renders  the 
two  kinds  of  graphic  outputs,  one  of  which  is  broadband  and  the  other  is  for  a  single 
frequency  data  file.  The  former  is  shown  in  Figure  3.16  (a)  and  the  latter  is  shown  in 
Figure  3.16  (b).  Note  that  if  multiple  frequencies  exist  in  the  data,  this  option 
automatically  creates  a  frequency  versus  depth  TL  field  plot. 


3  MMPE  2D  Graphic  Control  Frame 


MMPE  Output  File  Name      FreqPeoutbirv       _ 


Single  Radial   Single  Depth   Single  Interface  Source  Field  j 


Graphic  Output  of  Starting  Field 


Plot 


Help- 


n*  m 


Figure  3. 15  -  The  Option  "Source  Field"  in  Graphic  Control  Frame 
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Figure  3.16  -  Source  Field  Graphic  Output:  (a)  Broadband  Acoustic  Field;  (b)  Single  - 
Frequency  Acoustic  Field 
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IV.  OUTPUT  RESULTS 

A.  INTRODUCTION 

In  Chapters  II  and  III,  the  implementation  and  data  processing  of  the  MMPE  and 
WinMMPE  models  have  been  described.  In  this  chapter,  a  comparison  of  the  model 
output  data,  the  values  of  transmission  loss,  and  graphic  outputs  between  MMPE  and 
WinMMPE  are  made  in  order  to  confirm  the  accuracy  of  the  new  model.  The  efficiency 
of  WinMMPE  compared  with  MMPE  is  also  examined.  Specifically,  this  comparison  is 
based  upon  the  amount  of  memory  used  and  the  computational  run-times. 

As  examples  of  typical  cases,  two  sets  of  environmental  input  data  files  have  been 
defined.  The  first  case  is  based  on  a  simple,  range-independent  Pekeris  waveguide.  The 
water  column  is  isospeed  with  sound  speed  1500  m/s.  The  bottom  is  at  constant  depth 
250  m,  and  consists  of  a  homogeneous,  fluid  half  space  with  sound  speed  1600  m/s, 
density  1.2  g/cm3,  and  attenuation  0.15  dB/km/Hz.  The  source  is  at  a  depth  of  100  m  and 
transmits  a  CW  tone  of  100  Hz.  The  second  case  is  based  on  real  environmental 
measurements  taken  during  the  ONR-sponsored  Primer  experiment  off  the  east  coast  of 
New  Jersey  in  1996  (Lynch  et  al,  1997).  The  source  was  at  a  depth  of  270  m  and  the 
acoustic  energy  propagated  up  the  continental  slope  onto  the  shelf.  The  acoustic 
frequency  for  this  run  was  400  Hz.  Hereafter,  the  first  test  case  (Pekeris  waveguide)  will 
be  referred  to  as  Case  I,  and  the  second  test  case  (Primer)  will  be  referred  to  as  Case  EL 

As  shown  in  Table  4.1,  Case  I  has  only  one  profile  within  each  of  the 
environmental  input  data  files.  The  total  number  of  environmental  data  values  are  20.  By 
contrast,  Case  II  has  319,872  different  sound  speed  profiles  and  29,952  different 
bathymetry  data  points.  The  total  number  of  environmental  data  values  are  699664. 
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Test  Case 

Input  Data  Type 

Case  I 

Case  II 

Sound  Speed  Profile  Data  (2  data  values) 

1 

128radials  X 

49  ranges  X 

51  depths 

Bottom  Data  (2  data  values) 

1 

128  radials  X 
234  ranges 

Bottom  Property  Data  (7  data  values) 

1 

1 

Deep  Bottom  Data  (2  data  values) 

1 

1 

Deep  Bottom  Property  Data  (7  data  values) 

1 

1 

Table  4.1-  The  Number  of  Input  Data  Profiles 

Both  the  input  data  and  specification  of  parameters  within  "pefiles.inp"  decide  the 
matrix  of  the  output  data  from  the  model.  Specifically,  all  radials  defined  by  the 
environment  and  all  frequencies  defined  in  the  source  input  file  are  output.  However, 
only  the  number  of  depth  and  range  points  requested  between  certain  limits  defined 
within  "pefiles.inp"  are  stored.  Thus,  Case  I  has  a  two-dimensional  matrix,  as  shown  in 
Figure  4.1  (a),  of  size  103x200.  Case  n,  on  the  other  hand,  has  a  three-dimensional  data 
matrix,  as  shown  in  Figure  4.1  (b),  of  size  128x128x200. 


103  depthsteps 


(a) 


200  Range  Steps 


128Radia 


128  Depth  Steps 


(b) 


Figure  4. 1  -  Matrix  dimensions  of  output  data  for  (a)  Case  I  and  (b)  Case  II. 
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B.  COMPARISON  OF  OUTPUT  DATA 

The  Fortran  version  of  MMPE  primarily  uses  single  precision  floating  point 
values  with  each  data  occupying  4  bytes.  In  contrast,  the  Java  implementation  of 
WinMMPE  uses  8-byte,  double  precision  data  types.  There  are  two  reasons  for  using 
double  precision.  One  is  that  the  default  variable  type  in  Java  is  double  precision.  If  a 
developer  would  like  to  use  single  precision  variables  instead,  the  cast  operator  with  the 
"float"  command  must  be  added  in  the  definition  of  the  value.  Another  reason  to  use 
double  precision  is  that  more  precise  outputs  can  be  obtained.  However,  the  output  data 
from  both  models  is  single  precision. 

This  section  shows  the  comparison  of  output  data  from  the  MMPE  model  and  the 
WinMMPE  for  the  given  environmental  input  data.  Some  of  the  values  to  be  examined 
will  be  the  maximum  difference  of  the  PE  field  function  and  the  maximum  difference  of 
transmission  loss  (TL)  between  MMPE  and  WinMMPE  for  each  graphic.  In  addition, 
more  general  comparisons  of  the  graphic  outputs  from  each  model  are  made. 

1.  Single  Radial 

We  begin  by  examining  the  differences  between  single  radial  results.  The  PE 
field  function  as  defined  in  Eq.  (2.16)  is  a  complex  valued  vector.  As  shown  in  Table 
4.2,  the  maximum  amplitude  difference  of  the  PE  field  function  for  Case  I  is 
2.4780x10"*  while  the  maximum  difference  in  TL  is  0.0054  dB  re  lm.  These 
differences  are  negligible  and  can  be  attributed  to  the  differences  between  the  single 
precision  and  double  precision  implementations.  Figure  4.2  shows  the  "Single  Radial" 
graphic  outputs  from  MMPE  using  the  Matlab  post-processing  routines  and  WinMMPE. 

In  Case  n,  the  amplitude  of  the  maximum  difference  of  the  PE  field  function  is 
1.4950X10-4  while  the  maximum  difference  in  TL  is  7.5633  dB  re  lm  (Table  4.2).  This 
large  TL  difference  suggests  there  may  be  problems  between  the  solutions.  However, 
upon  closer  inspection,  it  is  found  that  all  of  the  significant  differences  occur  in  regions 
where  the  PE  field  function  is  extremely  small,  i.e.  at  discrete  locations  within  the  bottom 
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volume.  Table  4.3  lists  all  such  points  in  range  and  depth  where  the  TL  difference 
exceeds  3  dB.  In  total,  there  exist  8  points  over  128  depth  and  200  range  points.  Figure 
4.3  shows  the  "Single  Radial"  graphic  outputs  from  MMPE  and  WinMMPE  along  the 
central  radial.  These  plots  confirm  that  the  locations  of  the  discrepancies  are  in  regions 
of  extremely  low  field  strength.  Again,  this  could  be  due  simply  to  the  difference 
between  single  and  double  precision  implementations.  Overall,  these  plots  indicate  that 
WinMMPE  is  working  properly  and  producing  accurate  results. 


Test  Case 


Case  I 


Case  II 


Data 


PE 
Field  Function 


TL 

(dB  re  lm) 


PE 
Field  Function 


TL 
(dB  re  lm) 


Maximum  Difference 
( |MMPE  -  WinMMPE| ) 


2.4780X10" 


0.0054 


1.4950X10" 


7.5633 


Table  4.2  -  Maximum  difference  between  the  MMPE  Model  and  the  WinMMPE  Model 
for  Single  Radial. 


Table  4.3  -  F 


Depth  (m) 

Range  (km) 

TL  (dB  re  lm) 

219.5313 

14.3036 

3.0050 

341.4063 

18.9741 

7.5633 

363.2813 

20.4337 

3.7791 

391.4063 

20.4337 

4.4738 

307.0313 

23.9366 

5.3957 

297.6563 

42.6188 

3.0633 

269.5313 

46.4136 

4.3954 

257.0313 

49.3327 

3.4411 

Positions  where  T 

L  difference  exce 

:eds  3  dB  for  Case  II. 
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Figure  4.2  -  Single  Radial  graphic  output  for  Case  I:  (a)  from  MMPE;  (b)  from 
WinMMPE. 


47 


*  rtgi*?  Ho  I 


He    E.*     Toots    WM»o    Ht*> 


.  nix' 


20  30  40 

Rirge  [)  -h 


(a) 


&  lir-ii-  P«»J  (ivMLIM  4   .-Hv.+.Un} 


t^amsxissicw  loss  <:<J9  re  i  m;. 


,:i  »  411 

WHU    it",. 


Tl  ln#ofnwtton  tak*  »**>  ttv  Invj^J 


-J   Jul 


mo 


(b) 


Figure  4.3  -  Single  Radial  graphic  outputs  for  Case  II:  (a)  from  MMPE;  (b)  from 
WinMMPE. 
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2.  Single  Depth 

We  now  consider  the  data  produced  by  selecting  the  single  depth  output  option. 

The  maximum  amplitude  difference  of  the  PE  field  function  for  Case  I  is  2.1393X10"6 
while  the  maximum  difference  of  TL  is  0.0010  dB  re  lm.  These  values  are  recorded  in 
Table  4.4.  Again,  these  differences  may  be  considered  negligible.  Figure  4.4  shows  the 
"Single  Depth"  graphic  outputs  from  MMPE  and  WinMMPE  at  the  source  depth. 

For  Case  II,  the  maximum  amplitude  difference  of  the  PE  field  function  is  0.0016 
while  for  TL  the  difference  is  6.5112  dB  re  lm.  As  before,  the  positions  in  space  where 
the  TL  difference  exceeds  3  dB  are  determined  and  listed  in  Table  4.5.  In  this  case,  there 
exist  9  points  total  over  128  radial  and  200  range  points.  Figure  4.5  shows  the  "Single 
Depth"  graphic  outputs  from  MMPE  and  WinMMPE  at  the  source  depth.  As  expected, 
these  large  differences  only  occur  where  the  acoustic  field  has  extremely  low  values.  The 
similarity  of  the  plots  indicates  that  WinMMPE  is  providing  accurate  solutions. 


Test  Case 

Data 

Maximum  Difference 
( |MMPE  -  WinMMPE| ) 

Case  I 

PE 
Field  Function 

2.1393X10"6 

TL 
(dB  re  lm) 

0.0010 

Case  II 

PE 
Field  Function 

0.0016 

TL 
(dB  re  lm) 

6.5112 

Table  4.4  -  Maximum  difference  between  MMPE  and  WinMMPE  for  Single  Depth. 


Table  4.5  -  I 


Radial  (degree) 

Range  (km) 

TL  (dB  re  lm) 

-40.0781 

54.0033 

3.3533 

-23.9063 

57.5062 

3.3075 

-0.7031 

46.7055 

3.2034 

0 

46.4136 

3.0249 

11.2500 

51.6680 

5.6628 

16.1719 

27.4395 

4.2953 

26.0156 

11.3845 

6.5112 

30.2344 

25.6880 

5.0002 

40.0781 

14.8874 

4.9601 

Positions  where  T 

L  difference  exec 

jeds  3  dB  for  Case  II. 
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Figure  4.4  -  Single  Depth  graphic  outputs  for  Case  I:  (a)  from  MMPE;  (b)  from 
WinMMPE 
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Figure  4.5  -  Single  Depth  graphic  outputs  for  Case  II:  (a)  from  MMPE;  (b)  from 
WinMMPE. 
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3.  Single  Interface 

As  was  described  in  Chapter  HI,  there  are  two  radio  buttons  in  the  "Single 
Interface"  graphic  frame  which  allow  the  user  to  extract  data  either  from  the 
"Water/Bottom"  interface  or  the  "Bottom/Basement"  interface.  Because  the 
environments  examined  here  do  not  contain  the  deep  bottom/basement  interface,  we  will 
only  examine  the  "Water/Bottom"  interface  option  in  this  section.  Table  4.6  shows  that 
the  maximum  amplitude  difference  of  the  PE  field  function  for  Case  I  is  1.5451  xlO-6 
while  the  difference  in  TL  is  3.6884xl0~4dB  re  lm,  both  of  which  may  be  considered 
negligible.  Figure  4.6  shows  the  "Single  Interface"  graphic  outputs  from  MMPE  and 
WinMMPE  at  the  "Water/Bottom"  interface. 

For  Case  II,  the  maximum  amplitude  difference  of  the  PE  field  function  is  0.0017 
and  the  maximum  difference  in  TL  is  0.5484  dB  re  lm,  as  stated  in  Table  4.6.  Because 
we  are  explicitly  examining  values  along  the  interface,  we  do  not  encounter  the  extremely 
small  values  of  the  acoustic  field  that  were  problematic  before.  Furthermore,  the  largest 
differences  occur  in  the  TL  "nulls",  which  may  also  be  affected  by  the  double  precision 
implementation.  As  a  result,  even  these  differences  for  Case  II  may  be  considered 
negligible.  Figure  4.7  shows  the  "Single  Interface"  graphic  outputs  from  MMPE  and 
WinMMPE  at  the  "Water/Bottom"  interface. 


Test  Case 

Data 

Maximum  Difference 
( |MMPE  -  WinMMPEl ) 

Case  I 

Acoustic 
Field  Function 

1.5451X10^ 

TL 

(dB  re  lm) 

3.6884X10-4 

Case  II 

Acoustic 
Field  Function 

0.0017 

TL 

(dB  re  lm) 

0.5484 

Table  4.6  -  Maximum  difference  between  MMPE  and  WinMMPE  for  Water/Bottom 
Interface. 
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Figure  4.6  -  Water/Bottom  Interface  graphic  outputs  for  Case  I:  (a)  from  MMPE;  (b) 
from  WinMMPE. 
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Figure  4.7  -  Water/Bottom  Interface  graphic  outputs  for  Case  II:  (a)  from  MMPE;  (b) 
from  WinMMPE. 
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4.  Source  Field 

We  now  examine  the  differences  in  the  computed  starting  fields  at  zero  range. 
Table  4.7  shows  that  the  maximum  amplitude  difference  of  the  PE  field  function  for  Case 
I  is  2.5333  xlO-7  and  the  maximum  difference  in  TL  is  0.0664  dB  re  lm.  For  Case  II,  the 
maximum  amplitude  difference  of  the  PE  field  function  is  3.1593X10-6  while  the 
maximum  difference  in  TL  is  0.2548  dB  re  lm.  In  both  cases,  the  differences  may  be 
considered  negligible. 

Figure  4.8  displays  the  "Source  Field"  graphic  outputs  from  MMPE  and 
WinMMPE  at  the  initial  range  for  Case  I.  The  corresponding  graphic  outputs  for  Case  II 
are  provided  in  Figure  4.9. 


Test  Case 

Data 

Maximum  Difference 
( |MMPE  -  WinMMPE| ) 

Case  I 

Acoustic 
Field  Function 

2.5333X10"7 

TL 
(dB  re  lm) 

0.0664 

Case  II 

Acoustic 
Field  Function 

3.1593X10"6 

TL 
(dB  re  lm) 

0.2548 

Table  4.7  -  Maximum  difference  between  MMPE  and  WinMMPE  for  Source  Field. 
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Figure  4.8  -  Source  Field  graphic  outputs  for  Case  I:  (a)  from  MMPE:  (b)  from 
WinMMPE. 
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Figure  4.9  -  Source  Field  graphic  outputs  for  Case  II:  (a)  from  MMPE;  (b)  from 
WinMMPE. 
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C.  EFFICIENCY 

The  environment  defining  Case  I  was  very  simple  with  only  a  small  number  of 
input  parameters.  It  was  also  range-independent,  so  the  propagation  operator  only  needed 
to  be  defined  once  and  could  be  used  throughout  the  calculation.  Case  II,  on  the  other 
hand,  had  a  very  large  number  of  environmental  input  parameters  and  was  range- 
dependent.  This  required  an  update  of  the  propagation  operator  at  every  range  step. 
Therefore,  it  is  to  be  expected  that  Case  II  will  take  longer  to  run  than  Case  I. 
Furthermore,  the  amount  of  memory  used  for  Case  II  should  be  larger  than  for  Case  I. 

This  section  describes  the  efficiency  of  the  MMPE  Model  and  the  WinMMPE 
Model  as  both  a  propagation  model  and  a  graphic  tool  in  terms  of  main  memory  usage 
and  running  time.  All  testing  was  run  on  a  single  system  with  the  following 
specifications: 

•  Intel  Pentium  II  Processor,  366  MHz; 

•  Main  Memory,  392  Mbytes; 

•  Operating  System,  Microsoft  Windows  2000; 

•  Digital  Visual  Fortran  Compiler,  Standard  Edition,  Version  5.0.A; 

•  Matlab,  Version  5.3.1  (Rl  1.1)  with  Signal  Processing  Toolbox; 

•  Java2  SDK,  Standard  Edition,  Version  1.3.0; 

•  Forte  For  Java,  Community  Edition,  Version  1.0,  Update  Release  1. 

The  main  memory  on  this  system  is  adequate  to  avoid  the  need  for  swapping  with  the 
hard  disk.  Thus,  although  faster  systems  are  available  today,  this  should  provide  a 
reasonable  platform  to  compare  model  performance. 
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1.  Main  Memory 

As  mentioned  in  Chapter  I,  programmers  using  standard  Fortran  must  define  the 
size  of  data  arrays  at  the  beginning  of  the  program.  This  fixes  the  size  of  all  arrays 
throughout  the  calculation,  regardless  of  their  usage.  More  recent  versions  of  Fortran 
compilers  can  handle  dynamic  memory  allocation,  but  this  will  not  be  available  on  all 
platforms.  In  this  analysis,  therefore,  we  will  assume  the  Fortran  program  has  a  fixed 
size  of  data  arrays  for  input  data.  If  users  want  to  use  larger  input  data  sets,  the  program 
must  be  re-coded  to  accommodate  the  larger  data  array  sizes.  In  general,  when  a  very 
simple  data  set  is  used  which  only  occupies  a  fraction  of  the  fixed  allocated  space,  the 
rest  of  the  memory  is  wasted. 

In  Java,  however,  programmers  create  the  data  arrays  dynamically.  This  allows 
the  program  to  allocate  memory  space  only  when  it  needs  it  and  only  as  much  as  it  needs. 
This  provides  for  a  much  more  efficient  use  of  system  resources  than  the  Fortran  version. 
Even  though  the  Java  data  are  8-byte  double  precision  variables  as  compared  to  the 
mostly  4-byte  single  precision  variables  in  the  Fortran  version,  this  dynamic  allocation  of 
memory  is  expected  to  significantly  reduce  the  memory  requirements  of  the  program. 

The  amount  of  memory  used  in  each  version  of  the  model  is  shown  in  Table  4.8. 
These  values  were  obtained  from  the  "Task  Manager"  program  of  the  operating  system. 
Although  these  numbers  do  not  represent  exact  values,  they  do  provide  a  reasonable 
estimate  of  the  memory  usage  of  each  model.  For  MMPE,  both  Case  I  and  Case  II 
roughly  took  120  Mbytes.  The  memory  usage  of  WinMMPE  is  due  to  both  basic  model 
functions  (e.g.,  the  GUI  objects),  and  the  dynamic  allocation  for  run-time  variables.  The 
basic  functions  took  about  18  Mbytes.  The  dynamic  allocation  portion  took  about  4 
Mbytes  for  Case  I  and  17  Mbytes  for  Case  n. 
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Model 

Main  Memory 

Case  I 

Case  II 

MMPE  (Fortran) 

=  120  Mbytes 

=  120  Mbytes 

WinMMPE  (Java) 

=  1 8  Mbytes  (basic) 
+  4  Mbytes  (dynamic) 

=  18  Mbytes  (basic) 
+  17  Mbytes  (dynamic) 

Table  4.8  -  Main  memory  used  in  the  models  for  Case  I  and  Case  II. 

For  the  graphical  outputs  and  post-processing,  the  amount  of  memory  used  in 
WinMMPE  is  almost  the  same  as  that  used  in  Matlab.  This  is  because  Matlab  uses 
double  precision  and  dynamically  declares  the  size  of  data  arrays  with  given  sizes  from 
the  output  binary  file. 


2.  Running  Time 

Because  the  model  is  based  on  the  split-step  Fourier  algorithm,  a  significant 
fraction  of  the  run-time  is  due  to  the  fast  Fourier  transforms  (EFT)  between  depth  and 
vertical  wavenumber.  In  Case  I,  the  number  of  depth  points,  i.e.,  number  of  terms  in  the 
FFT,  is  1024.  In  Case  II,  the  number  of  depth  points  is  2048.  For  each  range  step,  two 
FFT's  must  be  performed.  Case  I  requires  222  range  steps  as  compared  to  15490  range 
steps  for  Case  II.  Furthermore,  for  multiple  radials  or  multiple  frequencies,  the  same 
number  of  FFT's  must  be  performed  for  each  range  calculation.  The  total  number  of 
FFT's  for  some  calculations  can  then  become  significant.  Table  4.9  shows  how  many 
times  the  FFT  module  in  the  model  is  called  for  each  case.  Note  that  Case  II  has  over 
8900  times  more  FFT  calls  than  Case  I. 


Tab 


Loop 


Radials 


Frequencies 


Ranges  x  2 


Total  =  Radials  x  Frequencies  x  Ranges  x  2 


Number  of  FFT  calls 


Case  I 


444 


444 


Case  II 


128 


30980 


3965440 


e  4.9  -  Number  of  FFT  calls  in  the  PE/SSF  algorithm  for  Case  I  and  Case  II. 
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Table  4.10  shows  the  run-times  of  the  models  for  Case  I  and  Case  II.  This 
measurement  was  achieved  by  watching  the  timer  in  the  operating  system.  In  some 
cases,  it  is  hard  to  tell  how  long  the  model  takes  to  complete,  especially  in  the  case  of 
small  input  data  and  relatively  few  range  steps.  For  Case  I,  the  MMPE  Model  spent  less 
than  a  second  and  the  WinMMPE  Model  spent  less  than  two  seconds.  Thus,  it  is  hard  to 
compare  the  true  relative  run-times.  Case  II,  however,  took  a  much  longer  time  to  run. 
The  MMPE  Model  took  about  3  hours  and  30  minutes.  In  contrast,  the  WinMMPE 
Model  took  6  hours  and  30  minutes.  It  appears,  therefore,  that  WinMMPE  takes  about 
twice  as  long  to  complete  a  calculation  as  MMPE.  This  factor  of  two  can  be  a  significant 
issue  for  very  large  calculations.  It  is  quite  possible  that  the  problem  is,  in  part,  due  to 
the  Java  code  method  of  calling  the  "getFFT(  )"  in  the  "FFT.java"  class  in  order  to 
implement  the  FFT.  The  Java  class  was  designed  like  the  Fortran  FFT  subroutine  in  the 
MMPE  Model,  which  is  based  on  a  standard  numerical  technique  (Press,  et  al.,  1992).  It 
is  recommended  that  the  FFT  routine  be  re-coded  in  an  object-oriented  programming 
style.  This  will  be  discussed  further  in  Chapter  V. 


Model 

Running  Time 

Case  I 

Case  II 

MMPE  (Fortran) 

<  lsec 

=  3hrs  30min 

WinMMPE  (Java) 

<  2  sec 

~  6hrs  30min 

Table  4. 10  -  Run-times  of  the  MMPE  and  WinMMPE  Models  for  Cases  I  and  II 


The  run-time  of  the  post-processing  and  graphic  presentation  is  listed  in  Table 
4.11.  For  Case  I,  the  Matlab  routines  spent  a  maximum  of  2  seconds  generating  the 
"Single  Radial"  output  while  the  WinMMPE  Model  spent  only  3  seconds.  The  run-times 
for  Case  II  are  nearly  the  same  for  the  two  models.  For  the  "Single  Depth"  graphic 
output,  the  run-time  of  the  WinMMPE  Model  is  typically  4  seconds.  The  Matlab 
routines,  on  the  other  hand,  consistently  took  about  41  seconds  to  complete.  For  the 
"Single  Interface"  graphic  output,  WinMMPE  generally  took  about  24  seconds  while 
Matlab  still  took  about  41  seconds  to  complete.  The  run-time  of  the  other  optional 
graphic  output,  "Source  Field,"  was  not  significant. 
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It  should  be  noted  that  these  run-times  do  not  account  for  the  time  taken  by  the 
user  to  run  the  propagation  model  and  the  post-processing  routines  separately.  Time  is 
also  taken  to  manually  enter  parameter  values  on  the  command  line  interfaces  for  the 
MMPE  Model  and  the  associated  Matlab  routines,  as  well  as  the  time  taken  to  fill  in  the 
parameter  fields  within  the  WinMMPE  Model.  The  WinMMPE  Model  at  least  has  the 
advantage  of  a  unified  program  within  which  GUI's  exist  to  control  the  propagation 
model,  post-processing,  and  graphic  presentations.  As  a  result,  the  operator's  skill  and 
familiarity  with  the  programs  can  affect  the  run-times. 


Graphic  Tool 

Run-Time 

Case  I 

Case  II 

Matlab 

Single  Radial:  <  2sec 

Single  Depth:  <  lsec 

Single  Interface:  <  lsec 

Source  Field:  <  lsec 

Single  Radial:  =  3sec 

Single  Depth:  =  4 lsec 

Single  Interface:  =  45sec 

Source  Field:  <  lsec 

WinMMPE  (Java) 

Single  Radial:  <  3 sec 

Single  Depth:  <  lsec 

Single  Interface:  <  lsec 

Source  Field:  <  lsec 

Single  Radial:  =  3sec 

Single  Depth:  =  4secs 

Single  Interface:  =  24secs 

Source  Field:  <  lsec 

Table  4.11-  Run-times  of  the  post-processing  and  graphic  presentations  for  Case  I  and 
Case  II 
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V.  CONCLUSIONS  &  RECOMMENDATIONS 

There  are  a  variety  of  underwater  acoustic  propagation  models  available  for  the 
general  user  these  days.  These  include  models  based  on  ray  theory,  normal  modes, 
wavenumber  integration  techniques,  and  parabolic  equation  models.  The  ONR- 
sponsored  web  site,  http://oalib.saic.com,  provides  copies  of  many  current  research 
models  used  in  the  professional  community.  However,  these  models  are  not  usually  very 
user-friendly  for  the  untrained.  None  of  the  models  examined  provide  a  convenient 
graphical  user  interface  and  few  have  been  developed  in  an  object-oriented  language.  To 
address  these  issues  and  provide  a  more  user-friendly  operating  environment,  a  Windows 
version  of  the  MMPE  Model,  referred  to  as  the  WinMMPE  Model,  was  developed. 

The  MMPE  Model  is  really  composed  of  two  parts.  The  propagation  component 
was  written  in  Fortran,  while  the  post-processing  and  graphical  presentations  were  treated 
using  Matlab  routines.  The  WinMMPE  Model,  on  the  other  hand,  incorporated  both 
parts  into  a  single  program  written  in  Java.  In  order  to  compare  the  accuracy  and 
efficiency  of  the  WinMMPE  Model  with  the  previous  MMPE  Model,  two  test  cases  were 
defined.  The  first  case  was  a  simple,  range-independent,  single  radial  Pekeris  waveguide 
with  very  few  environmental  input  parameters.  The  second  case  was  a  much  more 
complicated,  range-dependent,  3D  environment  with  multiple  radials  based  on  measured 
environmental  data  taken  near  the  Mid- Atlantic  Bight  off  the  coast  of  New  Jersey. 

A  comparison  of  the  data  indicated  that  WinMMPE  was  producing  results  nearly 
identical  to  MMPE.  The  largest  differences  were  observed  in  regions  where  the  acoustic 
field  strength  was  extremely  small  and  so  could  be  due  to  simple  round-off  error 
differences.  This  is  certainly  likely  since  the  default  data  type  in  the  Java  implementation 
of  WinMMPE  is  double  precision  while  the  default  type  in  the  Fortran  implementation  of 
MMPE  is  only  single  precision.  In  regions  where  the  acoustic  field  strength  was  large, 
the  differences  were  only  fractions  of  a  dB. 
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A  comparison  of  the  efficiency  of  each  model  in  terms  of  the  memory  usage  and 
run-time  was  also  presented.  Due  to  the  initialization  of  all  variable  array  sizes,  the 
Fortran-based  MMPE  Model  used  the  same  amount  of  memory  for  both  test  cases.  In 
contrast,  the  Java-based  WinMMPE  Model  dynamically  allocated  the  necessary  memory. 
The  MMPE  Model  always  required  about  120  Mbytes,  as  compared  to  the  roughly  25  - 
30  Mbytes  needed  to  run  the  propagation  component  of  the  WinMMPE  Model.  With 
respect  to  the  run-times,  is  seemed  that  the  WinMMPE  Model  generally  took  about  twice 
as  long  to  compute  the  propagation  solution  as  the  MMPE  Model.  This  could  be  due  in 
part  to  the  increase  from  single  precision  to  double  precision.  It  was  also  suggested  that 
the  FFT  routine  was  not  well  suited  for  implementation  within  Java.  Finally,  it  was 
found  that  there  was  essentially  no  difference  in  run-time  for  the  post-processing  and 
graphic  presentation  of  results. 

There  are  five  essential  recommendations  for  future  work.  The  first 
recommendation  concerns  the  FFT  routine.  As  suggested  above,  the  FFT  Java  class 
should  be  re-created  or  replaced  with  an  existing  FFT  Java  open  source  on  the  Internet. 
There  is  currently  no  FFT  Java  class  in  the  standard  library  provided  by  Sun 
Microsystems  Corp.  However,  more  and  more  class  libraries  are  becoming  available 
over      the      Internet.  One      of      them      is      provided      on      the      web      site 

http://www.fusion.kth.se/courses/jbone.  Complex  classes  must  be  combined  with  it  as 
we.  Therefore,  some  parts  of  the  algorithm  must  be  changed  since  the  treatment  of  real 
and  imaginary  data  must  be  combined  to  treat  a  complex  variable. 

Second,  there  are  additional  graphic  output  modules  that  need  to  be  added  in 

WinMMPE.   The  Matlab  routines  of  the  MMPE  Model  have  two  more  options  than  the 

WinMMPE  Model,  specifically  "Single  Range"  and  "Travel  Time".    The  first  of  these 

provides  another  view  of  the  data  to  help  the  user  evaluate  the  acoustical  environment. 

The  latter  option  is  of  great  interest  for  any  future  development  since  it  provides  for  the 

Fourier  synthesis  of  broadband  data  to  produce  predictions  of  pulse  propagation  in  the 

time  domain.     There  are  also  options  that  allow  the  user  to  produce  plane-wave 

beamformed  results,  thereby  providing  the  user  with  arrival  angle  versus  time  data. 
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Furthermore,  there  are  some  menu  items  in  the  graphic  frame  that  do  not  yet  work,  as 
mentioned  in  Chapter  III. 

Third,  the  problems  associated  with  the  "Input  File  Editor"  frame  were  described 
in  the  "Design  and  Functions"  section  of  Chapter  HI.  Even  though  that  frame  was  not  a 
major  work  in  this  thesis,  it  is  recommended  that  the  frame  be  changed  to  work  like  the 
well-functioned  standard  Windows  editor.  Additionally,  the  functions  of  the  menu  items 
in  the  "Edit"  menu  of  the  main  frame  should  be  upgraded  to  work  as  expected. 

Fourth,  it  would  be  a  great  service  to  the  general  acoustics  community  if  the 
WinMMPE  Model  was  available  over  the  Internet  using  Java  applets.  The  WinMMPE 
Model  might  provide  useful  information  about  the  underwater  acoustic  environment  for 
scientists  who  are  not  necessarily  acoustics  experts,  or  even  for  the  general  public 
interested  in  scientific  modeling  and  visualization.  It  could  even  be  used  by  schools  to 
interest  school  kids  in  science  and  the  ocean  environment. 

Finally,  3D  visualization  of  underwater  sound  propagation  can  be  an  interesting 
and  useful  project  in  the  future.  Holliday  (1998)  worked  on  real-time  3D  sonar  modeling 
and  visualization  with  a  ray-based  model.  His  work  provided  an  application- 
programming  interface  (API)  for  real-time  3D  sonar  modeling  and  visualization  of  sound 
propagation.  Future  work  could  consider  a  similar  3D  sonar  visualization  with  the 
WinMMPE  Model.  Figure  3.10  was  used  to  illustrate  the  scene  of  an  acoustic  field  for  a 
single  radial.  This  simple  example  was  not,  however,  a  true  3D  acoustic  field  since  it 
only  provided  a  single  plane  of  3D  structures  to  visualize.  Hopefully,  Figure  3.10  can  be 
a  starting  point  to  develop  the  3D  visualization  of  the  actual  3D  acoustic  data  provided  by 
WinMMPE. 
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